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Abstract 



Gluon jets are identified in e"'"e~ hadronic annihilation events by tagging two quark jets in 
the same hemisphere of an event. The gluon jet is defined inclusively as all the particles in 
the opposite hemisphere. Gluon jets defined in this manner have a close correspondence to 
gluon jets as they are defined for analytic calculations, and are almost independent of a jet 
finding algorithm. The charged particle multiplicity distribution of the gluon jets is presented, 
and is analyzed for its mean, dispersion, skew, and curtosis values, and for its factorial and 
cumulant moments. The results are compared to the analogous results found for a sample of 
light quark (uds) jets, also defined inclusively. We observe differences between the mean, skew 
and curtosis values of gluon and quark jets, but not between their dispersions. The cumulant 
moment results are compared to the predictions of QCD analytic calculations. A calculation 
which includes next-to-next-to-leading order corrections and energy conservation is observed to 
provide a much improved description of the data compared to a next-to-leading order calculation 
without energy conservation. There is agreement between the data and calculations for the 
ratios of the cumulant moments between gluon and quark jets. 
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1 Introduction 



Many experimental studies of quark jets have been performed at e^e~ colliders. Such studies 
are natural, since hadronic events in e"'"e~ annihilations above the T region and below the 
threshold for W^W~ production are believed to arise uniquely from the point-like creation of 
quark- ant iquark qq pairs. Production of the qq pair from a color-singlet point source allows 
the quark jets to be defined inclusively, by sums over the particles in an event or the event 
hemispheres. In contrast, conclusive experimental studies of gluon jets have been rare. This 
is because the creation of a gluon jet pair, gg, from a color singlet point source - allowing an 
inclusive definition analogous to that described above for quark jets - has been only rarely 
observed in nature.[] In most studies of gluon jets at e"'"e~ colliders, a jet finding algorithm 
is used to select an exclusive sample of three-jet qqg events. The same jet finder is used to 
artificially divide the particles of an event into a gluon jet part and two quark jet parts. In 
general, the results depend strongly on the algorithm chosen. Furthermore, use of a jet finder 
precludes a quantitative test of QCD analytic predictions for gluon and quark jet properties. 
For the analytic calculations, the gluon and quark jet characteristics are given by inclusive sums 
over the particles in color singlet gg and qq events, respectively, as described above. Thus, the 
theoretical results are not restricted to three-jet events defined by a jet finder and do not employ 
a jet finder to assign particles to the jets. 

In 0, a method was proposed for LEP experiments to identify gluon jets using an inclusive 
definition similar to that used for analytic calculations. The method is based on rare events of 
the type e+e~— > qqginci., in which the q and q are identified quark (or antiquark) jets which 
appear in the same hemisphere of an event. The object gmci., taken to be the gluon jet, is 
defined by the sum of all particles observed in the hemisphere opposite to that containing the 
q and q. In the limit that the q and q are collinear, the gluon jet gmci. is produced under the 
same conditions as gluon jets in gg events. The ginci. jets therefore correspond closely to single 
gluon jets in gg events, defined by dividing the gg events in half using the plane perpendicular 
to the principal event axis. First experimental results using this method were presented in 0. 

The results in Q were limited to the mean charged particle multiplicity values of gluon and 
quark jets. In this paper, we extend this study to include the full multiplicity distributions. 
The data were collected using the OPAL detector at LEP. For the quark jet sample, we select 
light quark (uds) event hemispheres, as in 0. Use of light quark events results in a better 
correspondence between the data and the massless quark assumption employed for analytic 
calculations, while use of event hemispheres to define the quark jets yields an inclusive definition 
analogous to that of the gluon jets ginci.- The multiplicity distributions of the gluon and quark 
jets are analyzed for their mean, dispersion, skew and curtosis values. In addition, we perform 
a factorial moment analysis of the gluon and quark jet multiplicity distributions in order to 
test the predictions of QCD analytic calculations |^ Q of those moments. 



^It is possible to identify a pure source of gg events in radiative T decays, such as in T(3S) — > 7x(j followed 
by Xb ^ gg 0; however, the jet energies are only about 5 GeV in this case, which limits their usefulness for jet 
studies. 



3 



2 Detector and data sample 



The OPAL detector is described in detail elsewhere 0. The present analysis is based on a 
sample of about 3 708 000 hadronic Z° decay events collected by OPAL from 1991 to 1995. 
Charged tracks measured in the OPAL central detector and clusters of energy measured in the 
electromagnetic calorimeter were selected for the analysis using the criteria given in [0. To 
minimize double counting of energy, clusters were used only if they were not associated with a 
charged track. Each accepted track and unassociated cluster was considered to be a particle. 
Tracks were assigned the pion mass. Clusters were assigned zero mass since they originate 
mostly from photons. To eliminate residual background and events in which a significant 
number of particles was lost near the beam direction, the number of accepted charged tracks 
was required to be at least five and the thrust axis of the event, calculated using the particles, 
was required to satisfy | cos(6'thrust)| < 0.9, where 6'thrust is the angle between the thrust and 
beam axes. The residual background from all sources was estimated to be less than 1%. 



3 Gluon jet selection 



For this study, a gluon jet is defined inclusively by the particles observed in an e"'"e~ event 
hemisphere opposite to a hemisphere containing an identified quark and antiquark jet, as stated 
in the introduction. The selection of inclusive gluon jets, gmci., is performed using the technique 
presented in [0]. More details concerning the motivation for the selection choices are given 
there. To select the ginci. gluon jets, each event is divided into hemispheres using the plane 
perpendicular to the thrust axis. Exactly two jets are reconstructed in each hemisphere, using 
the k± ("Durham") jet finder [^]. The results for the gluon jet properties are almost entirely 
insensitive to this choice of jet finder, as is discussed in and below in sections § and [77^ . 
Next, we attempt to reconstruct a displaced secondary vertex in each of the four jets. Displaced 
secondary vertices are associated with heavy quark decay, especially that of the b quark. At 
LEP, b quarks are produced almost exclusively at the electroweak vertex: thus a jet containing 
a b hadron is almost always a quark jet. To identify secondary vertices in jets, we employ the 
method given in fl^. Briefly, a secondary vertex is required to contain at least three tracks, 
at least two of which have a signed impact parameter value in the r-0 planeQ with respect to 
the primary event vertex, b, which satisfies b/ai, > 2.5, with af, the error of b. For jets with 
such a secondary vertex, the signed decay length, L, is calculated with respect to the primary 
vertex, along with its error, cr^. To be tagged as a quark jet, a jet is required to have a visible 
energy of at least 5 GeV and a successfully reconstructed secondary vertex with L < 2.0 cm 
and L/aL > 3.5.0 The visible energy of a jet is defined by the sum of the energy of the particles 
assigned to the jet. We refer to a hemisphere with two tagged jets as a tagged hemisphere. 

We next examine the angles that the two jets in a tagged hemisphere make with respect to 
the thrust axis and to each other. If the two jets are close together, or if one of the two jets is 
much more energetic than the other, it is very likely that one of the two jets is a gluon jet due 
to the strong kinematic similarity to an event with gluon radiation from a quark and because 

^Our coordinate system is defined so that z is the coordinate parallel to the e~ beam axis, r is the coordinate 
normal to the beam axis, (j) is the azimuthal angle around the beam axis and 6 is the polar angle with respect to z. 
^In 1^, we utilized a more stringent requirement of L/a^ > 5.0. 
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of the finite probability for a gluon jet to be identified as a b quark jet (see below). To reduce 
this background, we require the angle between the jets and thrust axis to exceed 10° and the 
angle between the two jets to exceed 50°. A last requirement is that the two jets lie no more 
than 70° from the thrust axis in order to eliminate jets near the hemisphere boundary. In total, 
324 events are selected for the final gluon jet ginci. sample.^ There are no events in which both 
hemispheres are tagged. 



We estimate the purity of this sample using the Jetset parton shower Monte Carlo ITT 



including detector simulation and the same analysis procedures as are applied to the data. 



The Jetset sample is a combination of events generated using version 7.3 of the program with 



the parameter values given in |jT3[ and of events generated using version 7.4 of the program with 
the parameter values given in 0. The initial Monte Carlo samples have about 3 000 000 events 
for version 7.3 and 4 000 000 events for version 7.4. The two Jetset versions yield results which 
are consistent with each other to within the statistical uncertainties and so we combine them. 
Using the Jetset events, the hadron level jets are examined to determine whether they are 
associated with an underlying quark or antiquark jet. To perform this association, the Monte 
Carlo events are also examined at the parton level. We determine the directions of the primary 
quark or antiquark from the TP decay after the parton shower evolution has terminated. The 
hadron jet closest to the direction of an evolved primary quark or antiquark is considered to 
be a quark jet. The distinct hadron jet closest to the evolved primary quark or antiquark not 
associated with this first hadron jet is considered to be the other quark jet. With the final 
cuts, Jetset predicts that both jets in the tagged hemisphere are quark jets with (80.5 ± 1.6)% 
probability, where the uncertainty is statistical: this is the estimated purity of the gmci. gluon 
jet sample. As an alternative method to estimate the gluon jet purity, we determine the fraction 
of events in the Monte Carlo ginci. sample for which the evolved primary quark and antiquark 
are both in the hemisphere opposite the ginci. jet: this yields the same estimate as given above. 
The Monte Carlo predicts that about 80% of the background is comprised of b events in which 
a gluon jet is mistakenly tagged as a b jet, while the other 20% is comprised about evenly of 
u, d, s and c events in which both a quark jet and a gluon jet are mistakenly tagged as b jets. 
The background events occur mostly when two tracks from the decay of a A or Kg hadron are 
combined with a third track to define a secondary vertex. 

Because we rely on displaced secondary vertices to identify quark jets, the ginci. jets in our 
study are contained in heavy quark events. The Monte Carlo with detector simulation predicts 
that about 95% of the events in the ginci. sample are b events. This reliance on b events is not 
expected to affect our results since the properties of hard, acollinear gluon jets do not depend 
on the event flavor according to QCD. More details are given in 0. 

We note that the ginci. jet tag rate, defined by the ratio of the number of ginci. jets to the 
number of events in the initial inclusive multihadronic event sample, is (8.74±0.49 (stat.)) x 10"^ 
for the data and (8.83 ± 0.36 (stat.)) x 10"^ for the Monte Carlo. Thus the Monte Carlo 
reproduces the measured tag rate well. 

The mean energy of the gluon jets, {E)^^^^^ , is less than the beam energy because the 
two quark jets against which ginci. recoils are not entirely coUinear. The mean visible energy 



^Because we use a somewhat different event selection here, and because our data have been reprocessed using 
improved detector cahbrations since the time of our previous pubhcation only 60% of these events are in 
common with the ginci. sample in H. 
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of the ginci. jets, corrected for the effects of the detector and initial-state photon radiation, is 
41.8±0.6 (stat.) GeV. As an alternative method to estimate the gluon jet energy, we employ the 
technique of calculated jet energies for massive jets. A jet direction is determined for the gluon 
jet by summing the momenta of the particles in the gmci. hemisphere. The angles between the 
ginci. jet and the two jets in the tagged hemisphere are used in conjunction with the measured jet 
velocities to calculate the jet energies assuming energy-momentum conservation.^ The velocity 
of a jet is given by the magnitude of its visible 3-momentum divided by its visible energy. Using 
this method, the mean gluon jet energy is determined to be 41.5 ± 0.3 (stat.) GeV, which is 
consistent with the result given above for the visible energy. In this paper, we choose to use 
the visible ginci. energy, rather than the calculated energy, because it corresponds more closely 
to the jet energy as it is defined for the uds hemisphere quark jets. The difference between the 
mean visible and calculated jet energies is used to define a systematic uncertainty. The mean 
energy of the gluon jets in our study is therefore {E)g.^^^ =41.8 ± 0.6 (stat.) ± 0.3 (syst.) GeV. 



4 Monte Carlo comparison of gmci. cind gg jets 

Our analysis of gluon jets is based on the premise that ginci. jets from e+e^ annihilations are 
equivalent to hemispheres of gg events produced from a color singlet point source, with the 
hemispheres defined by the plane perpendicular to the thrust axis. Although high energy gg 
events are not available experimentally, they may be generated using a QCD Monte Carlo event 
generator. The viability of our premise can be tested by comparing the Monte Carlo predictions 
for gg event hemispheres and ginci. jets. Such a comparison has already been presented in |^ 
for the mean charged particle multiplicity values, (ndi.) (see also [0)- Here, we extend this 
comparison to the full multiphcity distribution, P(nch.) versus rich., with P(nch.) the probability 
that an event will be observed with a charged particle multiplicity rich.- 

The points with error bars in Figure |l| show the prediction of the Herwig parton shower 
Monte Carlo |]T5|, version 5.9, for the charged particle multiplicity distribution of ginci. jets. 
The uncertainties are statistical. The parameter set we use is the same as that given in 0] 
for Herwig, version 5.8, except that the value of the cluster mass cutoff CLMAX has been 
increased from 3.40 GeV/c^ to 3.75 GeV/c^ to improve the model's description of (ndi.) in 
inclusive hadronic Z° decays. The e"'"e~— > qq ginci. events were generated using a center-of-mass 
(cm.) energy, -Ec.m., of 91.2 GeV to correspond to the data. The ginci. identification was 
performed using the same procedure as is described for the data in section |^, except that the 
two quark jets against which the ginci. jet recoils were identified using the Monte Carlo method 
described in section ^ In particular, the angular cuts on the directions of the quark jets with 
respect to the thrust axis and to each other have been applied. The resulting mean energy of 
the Monte Carlo ginci. jets is 41.2 GeV with a negligible statistical uncertainty. 

Shown by the solid histogram in Figure |I| is the prediction of Herwig for gg event hemi- 
spheres. The gg events were generated using a cm. energy of 82.4 GeV so that the hemisphere 
energies are the same as for the ginci. jets. It is seen that the results for the ginci. jets and the gg 

^This method results in a better estimate of the ginci. jet energy than the method assuming massless jets 
which we employed in [^i the Monte Carlo without detector simulation yields identical results for the visible 
and calculated ginci. jet energies if the massive formula is used, whereas the calculated energy is about 2 GeV 
smaller than the visible one if the massless formula is used. 
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event hemispheres are essentially identical. This establishes the viability of our method, con- 
firming the results of [0]. Similar agreement between the predicted multiplicity distributions of 
ginci. jets and gg event hemispheres is obtained if Jetset is used to generate the samples rather 
than Herwig, or if the JADE-EO [0 or cone [0 jet finder is used to identify the quark jets for 
the ginci. jet selection, rather than the k± jet finder. 



5 uds quark jet selection 

The uds quark jets in our study are defined inclusively, by summing the particles observed in 
an event hemisphere opposite to a hemisphere containing an identified uds jet. Since there are 
only 324 gluon jets in our study, it is not necessary to use the entire data sample of about 
3 708 000 events mentioned in section || for the uds jet analysis. Instead, we base the uds jet 
selection on an initial sample of about 396 000 hadronic annihilation events with cm. energies 
within 100 MeV of the Z° peak. 

To select the uds jets, we divide each event into hemispheres using the plane perpendicular 
to the thrust axis. Selection criteria are applied to each hemisphere separately using charged 
tracks that appear in a cone of half angle 40° around the thrust axis. The reason for the restric- 
tion to tracks which lie within 40° of the thrust axis is to avoid using tracks near the hemisphere 
boundary. An algorithm is applied to identify charged tracks which are consistent with arising 
from photon conversions ]TB|]. Removing such tracks from consideration, the number of tracks 
in the cone which have a signed impact parameter significance, b/ab, greater than 1.5 is de- 
termined. A hemisphere is tagged as containing a uds jet if the number of tracks in the cone 
which have b/ab > 1-5 is zero. In total, 188 288 hemispheres are tagged. This number includes 
30 303 events for which both hemispheres are tagged. The estimated uds purity of this sample, 
obtained by treating Jetset events with detector simulation in the same manner as the data, 
is 80.9%, with a negligible statistical uncertainty. The Monte Carlo predicts that about 70% 
of the background events are c events and that 30% are b events. The uds jet tag rate, defined 
by the ratio of the number of tagged uds jets to the number of events in the initial inclusive 
multihadronic event sample, is 0.399 ±0.001 (stat.) for the data and 0.418 ±0.001 (stat.) for the 
Monte Carlo. The difference of about 2% between the uds jet tag rates of data and Monte Carlo 
implies a small deficiency in the simulation of the event characteristics, which is accounted for 
in our evaluation of systematic uncertainties (section [77^ ). The corrected energy of the uds jets 
is given by the beam energy, 45.6 GeV, with essentially no uncertainty. 

For purposes of comparison. Figure |I| includes the prediction of Herwig for uds qq event 
hemispheres, generated using the same cm. energy that is used to generate the gg event sample. 



6 Corrections 



The measured charged particle multiplicity distributions of the gmci. and uds jets are corrected 
in two steps, following the method presented in |jl9[. In the first step, the data are corrected for 
experimental acceptance, resolution, and secondary electromagnetic and hadronic interactions 
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using an unfolding matrix [IS]. This matrix is constructed using Jetset events, including full 
detector simulation and the same selection criteria as are applied to the data. In the second step, 
the data are corrected for event acceptance and the effects of initial-state photon radiation using 
bin-by-bin multiplicative factors. The bin-by-bin corrections are derived using two different 
Jetset samples. The first sample, based on inclusive TP hadronic decays, includes initial-state 
photon radiation and the same event acceptance criteria as the data, but not detector simulation 
(they are the generator level input to events which have been processed through the detector 
simulation and which have been selected using the same selection criteria as are applied to the 
data). The second sample does not include initial-state photon radiation, event acceptance, 
or detector simulation and treats all charged and neutral particles with mean lifetimes greater 
than 3 x 10^^° s as stable: hence charged particles from the decays of Kg and weakly decaying 
hyperons are included in the definition of multiplicity. For the correction of the gluon jet 
data, inclusive TP events are used for the second sample. The quark jets in this sample are 
identified with Monte Carlo information using the method discussed in section ^: otherwise 
the ginci. sample is obtained in the same manner as is described in section |^ for the data. For 
the correction of the uds jet data, the jets of the second sample are defined by the particles in 
the hemispheres of uds events. The multiplicative correction factors are obtained by taking the 
ratios of the predictions from the second sample to those from the first one. Thus, the bin-by- 
bin corrections account not only for detector response and initial-state radiation but also for 
the background to the gmci. and uds jet data. The corrections applied to the data are generally 
moderate or small. For example, Jetset predicts the mean multiplicity value (section |7. 1| ) of 
ginci. jets to be only 6% larger at the generator level than it is at the level which includes 
detector simulation and the experimental selection criteria. The corresponding difference for 
uds jets is —4%. 



7 Results 



The corrected charged particle multiplicity distributions are presented in Figure Numerical 
values for these data are listed in Table 0. Shown in comparison to the data are the generator 
level predictions of Herwig 5.9 and Jetset 7.4. The Monte Carlo results for the ginci. jets are 
obtained in the manner described in section ^. The two models are seen to provide a generally 
adequate description of the measurements except that the uds jet distribution predicted by 
Herwig is shifted towards lower values of multiplicity than are observed experimentally (Fig- 
ure 11(b)). Statistical uncertainties were estimated for the ginci. results using 100 Monte Carlo 
samples of ginci. jets at the generator level, each with approximately the same event statistics 
as the data sample. The statistical uncertainty for each result (e.g. a multiplicity bin in Fig- 
ure 11(a) or a factorial moment measurement, see section |7!^ ) was set equal to the RMS value 
found for the 100 samples. The same method was used to evaluate statistical uncertainties for 
the uds jets. The matrix corrections introduce correlations between the bins of the corrected 
multiplicity distributions. The correlations are generally strong between a bin and its nearest 
one or two neighbors on either side but can extend with smaller strength to four or five bins 
away. These correlations smooth out bin-to-bin statistical fluctuations. This effect is partic- 
ularly noticeable for the gluon jet distribution (Figure |^(a)) because of the relatively small 
number of events in the ginci. jet sample. 
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7.1 Mean, dispersion, skew and curtosis values 



We determine the mean {rich), dispersion D = ^ {n\) — {uch.y, skew 7 = {{n ch. — {n ch.)Y) / 
and curtosis c = [{{{rich. — {in ch)Y) / D'^) — 3] values of 41.8 GeV ginci. gluon jet and 45.6 GeV 
uds quark jet hemispheres to be: 

(1) 
(2) 
(3) 
(4) 



where the first uncertainty is statistical and the second is systematic. These results are shown 
in Figures ^ and ||. The systematic uncertainties are discussed below in section [773[ The results 
for (nch.)ginei (^ch.) uds hemis. are consistent with those presented in pO| . 





= 14.32 ±0.23 ±0.40 


ch. ) uds hemis. 


= 10.10 ±0.01 ±0.18 


6 incl. 


= 4.37 ±0.19 ±0.26 


Duds hemis. 


= 4.298 ±0.008 ±0.098 


7gincl. 


= 0.38 ±0.13 ±0.18 


Tuds hemis. 


= 0.822 ±0.007 ±0.044 


^ § incl. 


= 0.18 ±0.34 ±0.30 


C uds hemis. 


= 0.98 ±0.03 ±0.11 , 



Figures | and | include the Herwig and Jetset predictions for ginci. jets, shown by the cross 
and diamond symbols, respectively. Also shown, by the finely-dashed and solid horizontal lines, 
are the Monte Carlo predictions for gg event hemispheres, generated to have the same energy 
as the ginci. jets. It is seen that the results for the gg hemispheres and ginci. samples agree 
well for both Herwig and Jetset (compare the cross symbols to the finely-dashed lines and the 
diamond symbols to the solid lines), which is consistent with Figure |I| and the discussion in 
section |[ The Monte Carlo predictions for gluon jet properties are seen to agree well with the 
data, except that the Jetset prediction of the mean multiplicity (?T-ch.)ginci somewhat above 
the measured value (Figure ^a)). 

Also shown in Figures |^ and ^, by the coarsely-dashed and dash-dotted horizontal lines, are 
the predictions of Herwig and Jetset for uds event hemispheres. These predictions are shown 
for two different values of jet energy: £'jet=45.6 GeV (i?c.m.=91.2 GeV), corresponding to the 
energy of the measured uds jets, and i?jet=41.8 GeV (i?c.m.=83.6 GeV), corresponding to the 
energy of the ginci. jets. The steps in Figures]^ and || between the predictions for i?c.m.=91.2 GeV 
and Ec.m.=83.6 GeV uds hemispheres therefore indicate the Monte Carlo corrections for quark 
jets to account for the difference in energy between the uds and ginci. samples. Comparing the 
Monte Carlo predictions for i?c.m.=91.2 GeV to the data, it is seen that Jetset provides a good 
overall description of the uds jet properties. Herwig's predictions for the uds jet properties 
are also in reasonable agreement with the data, with the exception of the mean multiplicity 
('^ch.) uds hemis. (Figure |](a)), which is somewhat too low as was already noted in connection with 
Figure @(b). 

We also determine the ratios between the gluon and quark jet results since common system- 
atic uncertainties will partially cancel. Before forming these ratios, it is necessary to account 
for the different energies of the two samples: the gluon jets have a mean energy of 41.8 GeV 
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while the uds jets have a mean energy of 45.6 GeV. To correct the quark jet (^ch.) value for 
this difference in energy, we follow the method in ||^ and employ the QCD analytic formula 
for the evolution of the mean event multiplicity in e+e~ annihilations [^. This QCD result is 
known to describe the energy evolution of the mean charged particle multiplicity in inclusive 
e"'"e~ annihilation events with good accuracy [^. Assuming the number of active quark flavors, 
nj, to be five, the QCD evolution formula predicts the mean multiplicity of 41.8 GeV quark 
jet hemispheres to be (3.6 ± 0.2)% smaller than for 45.6 GeV quark jet hemispheres, where 
the uncertainty results from the maximum variation found by using the jet energies (41.8 GeV 
and 45.6 GeV) rather than the event energies (83.6 GeV and 91.2 GeV), nf=3 rather than 
nf=5, and varying the value of within its allowed range ||2^.Q Virtually the same result 
is obtained if the evolution formula is evaluated using the fitted values given in [|3 for the 
strong coupling strength and the overall normalization. Applying a multiplicative correction of 
0.964 to the result presented above for (^ch.) uds hemis. yields ('^ch.)uds\^mis.=9.74 ± 0.01 (stat.). 
Our result for the multiplicity ratio Vch. between 41.8 GeV gluon and quark jets is therefore: 

rch. = ^^wIIgcV = l-^''! ^ 0-0^4 (stat.) ± 0.043 (syst.) . (5) 

da.) uds homis. 

This result is consistent with our previous result but has substantially reduced uncertain- 
ties. Furthermore, the analytic prediction in p5[ is in general agreement with this measurement. 



For purposes of comparison, the predictions of Herwig and Jetset are rch. =1-537 ± 0.002 and 
1.539 ± 0.002, respectively, where the uncertainties are statistical. 

For the dispersion, skew and curtosis values, we account for the difference in energy between 
the gluon and quark jet measurements using the Monte Carlo predictions. The Monte Carlo is 
known to provide a good description of the energy evolution of the dispersion of the multiplicity 



distribution in inclusive e"'"e~ hadronic events (e.g. see [|T9|, |2^, |J]), making it plausible that its 
predictions for the energy evolution of skew and curtosis are also reliable. Jetset predicts the 
dispersion D to be (3.6 ± 0.2)% smaller for 41.8 GeV uds hemispheres than for 45.6 GeV uds 
hemispheres (Figure ^(b)). The corresponding results for skew and curtosis are (2.5 ± 0.5)% 
and (4.5 ± 2.2)%, respectively (Figure §(a) and (b)). The uncertainties for these values are 
given by the maximum variation found by using Herwig rather than Jetset, nf=5 rather than 
n/=3, and varying the value of Alla by its uncertainty 0. Applying corrections of 0.964, 0.975 
and 0.955 to the uds jet dispersion, skew and curtosis measurements given above, respectively, 
yields Z^^kSs. = 4.143 ± 0.008 (stat.), 7u'ds'h?^s. = 0.802 ± 0.007 (stat.), and c^^fS, = 
0.933 ± 0.027 (stat.). The ratios between the dispersion, skew and curtosis values of 41.8 GeV 
gluon and quark jets are therefore: 

= J^iTceV = 1-055 ± 0.046 (stat.) ± 0.055 (syst.) (6) 



n 41.8 GeV 
uds hemis. 



7, 



S incl. 



^41.8 GeV 
/ uds hemis. 

^ S incl. 

„41.8GeV 
uds hemis. 



0.47 ±0.16 (stat.) ±0.21 (syst.) (7) 
0.19 ±0.37 (stat.) ±0.33 (syst.) . (8) 



® Jetset and Herwig predict reductions of 3.5% and 3.6%, respectively, in the value of (rich.) for uds hemi- 
spheres with i?c.m.=83.6 GeV compared to those with Ec,i^=91.2 GeV: thus the Monte Carlo predictions for 
the energy correction are essentially the same as that obtained from the QCD evolution formula. 

''There is a shift of 0.081 between the central values of our previous result (rch. =1-552 ± 0.073 (stat.+syst.)) 
and our current result; this shift is due primarily to our reevaluation of the gmci. energy value (section |^) and 
thus to a reevaluation of the correction to account for the difference in energy between the uds and ginci. jets. 
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The dispersions of the gluon and quark jet multiphcity distributions are therefore almost equal, 
despite the large difference between their mean values (relation (^). The quark jet distribution 
is more skewed, i.e. asymmetric, than the corresponding distribution for gluon jets. Quark jets 
are observed to have a larger curtosis value than gluon jets, implying that their distribution is 
more non-gaussian in shape (the peak is higher and the tails are broader than a gaussian with 
the same mean and dispersion) than is the case for gluon jets. We note that the deviations 
of the ratios and from unity are only 2.0 and 1.6 standard deviations of their total 
uncertainties, however. Herwig predicts r£)=1.086 it 0.003, r^=0.60 it 0.01 and rc=0.30 ± 0.03, 
where the uncertainties are statistical. The corresponding results from Jetset are 1.110 ±0.003, 
0.55 ±0.01 and 0.31 ± 0.03. 



7.2 Factorial, cumulant and Hq moments 

Factorial moments provide a standard means to characterize the fluctuations of a distribution 
about its mean value (cf. ^, ^). Factorial moments are less subject to bias from random 
statistical fluctuations than "ordinary" central moments, as is discussed in p6|. Various QCD 
analytic calculations have been performed for the factorial moments of the multiplicity distribu- 
tions of separated gluon and quark jets |^. A factorial moment analysis of our data permits 
a test of these QCD calculations for the first time (see section |^). Currently, QCD predic- 
tions do not exist for the dispersion, skew and curtosis values presented in section [7.1| . Certain 
combinations of factorial moments are directly related to these three quantities, however, as is 
discussed below. 



The normalized factorial moment of rank g, Fg, is defined by p6| , p7| : 

_ (n(n- l)---(n-g± 1)) 

= ' 

where q > 1. An equivalent characterization is given by the cumulant factorial moments. 
Kg , which can be obtained from the Fg moments recursively: 



(a -IV 

K,^Fg-J2—^ F„ , (10) 

with the condition Ki=l. The Kg moments have been shown to be more sensitive to de- 
tailed features of the multiplicity distribution than the Fg moments they are more sensi- 
tive to higher order QCD corrections and to differences between QCD and phenomenological 
parametrizations such as the negative binomial distribution. Besides factorial and cumulant 
moments, it has become standard to consider the ratio of cumulant to factorial moments. 



denoted Hg, which appear naturally in the solution of the QCD equations 



H,^^ . (11) 

An analysis of the Hg factorial moments of the inclusive multiplicity distribution in multi- 
hadronic Z" decays has recently been presented in pQl BO]. 



In Tables ^ and ^, we present our measurements of the Fg, Kg and Hg factorial moments of 
ginci. and uds jets. Our results are given for ranks 2 < g < 5. {Fi=Ki=Hi=l trivially.) These 
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ranks correspond to those for which theoretical predictions have been pubhshed The 
results for Fg and Kq are shown in Figures ^ and The results for Hq are qualitatively similar 
to those shown in Figure |^ for Kq and so are not shown in addition. Figures ^ and ^ include 
the predictions of Herwig and Jetset. The Monte Carlo information is presented in the same 
manner as in Figures ^ and |^. It is seen that the Monte Carlo results for the ginci. and gg event 
hemispheres agree well with each other, i.e. the cross symbols agree with the finely-dashed lines 
and the diamond symbols agree with the solid lines to within differences that are consistent 
with statistical fluctuations. It is also seen that the Monte Carlo corrections to uds jets to 
account for the difference in energy between the ginci. and uds samples (the steps in the center 
of the coarsely-dashed and dash-dotted curves in Figures |^ and |^) are moderate in comparison 
to the experimental uncertainties. 

To gain insight concerning the physical interpretation of the factorial and cumulant mo- 
ments, and to help relate the measurements shown in Figures ^ and ^ to those shown in 
Figures |^ and ^, we generated 100 samples of ginci. jets using the Jetset Monte Carlo, each 
with event statistics similar to that of the data. We used these 100 generator level samples to 
calculate the correlation coefficients between the Fq moments, the Kq moments and the mean, 
dispersion, skew and curtosis values. The resulting correlation matrix is presented in Table |^. 
Similar results were found in an analogous study of uds quark jets. We also selected 24 indepen- 
dent samples of uds events which included detector simulation and the same selection criteria as 
the data. We processed these 24 samples using the correction procedure described in section |] 
and determined the correlations between the corrected results. The detector level study was 
not repeated for ginci. jets because of inadequate Monte Carlo event statistics. The resulting 
correlation matrix (corresponding to Table |^) was found to be very similar to that obtained 
using the 100 generator level samples, from which we conclude that the correction procedure 
does not introduce significant correlations between the variables. To facilitate the discussion 
in the next paragraph, boxes have been drawn around correlation coefficients in Table |^ which 
have magnitudes of 0.90 or larger. 

From Table ^, it is seen that there is a high degree of statistical correlation between factorial 
moments of different rank: the correlation coefficients between the various Fq moments lie 
between 0.92 and 0.99. In contrast, only modest or small correlations are observed between the 
cumulant moments. The largest correlation coefficient in this case, between the Kq moments 
with ranks g=3 and (3'=4, is only 0.43. Thus, the results shown for different ranks in Figure ^ are 
largely independent of each other, while those shown in Figure § are not. Table ^ also establishes 
that the cumulant moments of ranks 2, 3 and 4 are strongly correlated with dispersion, skew 
and curtosis, respectively (correlation coefficients of 0.97, 0.99 and 0.90). Furthermore, with the 
exception of the correlation between K^^ and curtosis (coefficient of 0.78), the other correlations 
of the cumulant moments with dispersion, skew and curtosis are moderate or small. Therefore, 
the cumulant moments of ranks 2, 3 and 4 are directly related to dispersion, skew and curtosis, 
respectively. Algebraically, the relationships are: 



^2 = - 1—, (12) 




^^-3 ^ .(7^1 - x^It^I + (13) 
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Thus, the i^'g moments with q=2, 3 and 4 are essentially equivalent to dispersion, skew and 
curtosis, but are in a form for which QCD analytic calculations have been presented (see 
section ||). In contrast, the Fq moments exhibit a strong correlation with dispersion but not 
with skew or curtosis, as is seen from Table |^. 

In Table |, we present measurements of the ratios rp^, tk^ and between the Fg, Kg and 
Hq factorial moments of 41.8 GeV gluon and quark jets. To obtain these results, the quark jet 
values in Table S were corrected for the difference in energy between the uds and ginci. samples 



using the method described in section |7]1] for the dispersion, skew and curtosis, i.e. using the 
Jetset predictions (e.g. Figures |]andP), with a systematic uncertainty evaluated as is described 
in section |7]l|. These corrections typically lie between 0.95 and 0.99. The ratios rp^, tk^ and 
THq are then formed by dividing the factorial moments of gluon jets (Table ^) by these corrected 
quark jet results. Our measurements of the cumulant moment ratios r^q are shown in Figure |^. 
It is seen that the gluon and quark jet cumulant moments differ by about a factor of three for 
q=2 and by an even larger amount for the higher moments. From relations ([T^)-(|l^, it is seen 
that part of this difference can be attributed to the difference between the mean values (nch.) 
of gluon and quark jets (relation (||)). The results are well reproduced by the predictions of 
Herwig and Jetset, shown by the dashed and solid horizontal lines in Figure |^. 



7.3 Systematic Uncertainties 

To evaluate systematic uncertainties, the analysis was repeated with the following changes 
relative to the standard analysis. There were no significant changes in the number of selected 
events or in their estimated purities compared to the standard results unless otherwise noted. 



1. Charged tracks alone were used for the data and for the Monte Carlo samples which 
include detector simulation, rather than charged tracks plus unassociated electromagnetic 
clusters. 

2. Herwig was used to determine the correction matrix and bin-by-bin correction factors, 
rather than Jetset. 

3. The particle selection was varied, first by restricting charged tracks and electromagnetic 
clusters to the central region of the detector, | cos(6')| < 0.70, rather than | cos(6')| < 0.94 
for the charged tracks and | cos(^)| < 0.98 for the clusters, and second by increasing the 
minimum momentum of charged tracks, p^^^ , from 0.10 GeV/c to 0.20 GeV/c. 



4. The gluon jet selection was performed using the JADE-EO [|T^] and cone |T^ jet finders to 
define the tagged quark jets, rather than the /c_l jet finder: 320 and 246 gmci. jets resulted, 
respectively, of which 88% and 76% were in common with the events of the standard ginci. 
sample. 

5. The geometric conditions for the gluon jet selection were varied, first by requiring the 
angle between the two jets in the tagged hemisphere to exceed 65°, rather than 50°, and 
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second by requiring the two tagged quark jets to lie within 65° of the thrust axis, rather 
than 70°. 

6. At least one track with a signed impact parameter significance greater than 2.5 was 
required to be present in the displaced secondary vertices used to tag quark jets for the 
ginci. identification, rather than at least two; the ginci. sample increased to 1127 jets, while 
its estimated gluon jet purity decreased to 55.8%. 

7. The gluon jet sample was restricted to events collected within 100 MeV of the Z° peak. 

8. uds jets were tagged using charged tracks that appeared within a cone of half angle 70° 
around the thrust axis, rather than 40°. 

9. The maximum signed impact parameter significance of tracks used for the identification 
of uds jets was increased from 1.5 to 2.5. 

10. For the ratios of mean multiplicity, dispersion, skew and curtosis Vch., rn, and Tc, and 
for the ratios of factorial moments rp^, r^q and Vh^, the energy to which the quark jet 
results were corrected was varied by the total uncertainty of the ginci. jet energy (section 
also, for these same quantities, the correction factors to account for the difference between 
the uds and gmci. jet energies were varied by their uncertainties (sections ^[l] and [7^ ). 



The differences between the standard results and those found using each of these conditions 
were used to define symmetric systematic uncertainties. For items 3, 4, 5 and 10, the larger of 
the two described differences with respect to the standard result was assigned as the systematic 
uncertainty. For item 2, the difference with respect to the standard result was multiplied by 
2/ a/T2 [^] since Herwig represents an extreme choice of hadronization model compared to Jet- 
set. For the uds jet differential multiplicity distribution (Figure ||(b) and Table |l]), we evaluated 
the systematic terms involving p^'j^ (item 3 in the above list) using the procedure described 
in |T9[ : the corrected distribution was parameterized using polynomials, the parameterized dis- 



tribution was shifted along the multiplicity axis so that its mean coincided with the mean of 
the standard result, and the systematic uncertainty was defined bin- by-bin by the difference 
between the shifted and standard distributions. 



The uncertainties were added in quadrature to define the total systematic uncertainty. For 
the differential multiplicity distributions (Figure || and Table [^), the systematic uncertainty 
evaluated for each bin was averaged with the results from its two neighbors to reduce the 
effect of bin-to-bin fluctuations (the single neighbor was used for bins on the endpoints of the 
distributions). The largest systematic terms for the ginci. jet measurements were generally found 
to arise about equally from items 3, 4 and 5 in the above list. The largest systematic terms 
for the uds jet measurements were generally found to arise from item 3 and, to a lesser extent, 
from item 1. For the ratios of the gluon to quark jet measurements, the largest systematic 
terms generally arose from items 1 and 3-5. As an illustration. Table ^ provides a breakdown 
of the systematic uncertainties evaluated for (nch.)ginci > ('^ch.)udshemis., and Vch.- 

From Table ^, it is seen that using the JADE-EO or cone jet flnders to identify the tagged 
quark jets for the ginci. sample (item 4 in the above list), yields rch.=l-488 ± 0.024 (stat.) or 
rch.=l-448 ± 0.028 (stat.), respectively, which differ by less than 2% from the standard result 
of Teh. =1-471 ± 0.024 (stat.). This emphasizes that our results are almost independent of the 
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choice of the jet finding algorithm. In contrast, results for Vch. based on exclusive samples 
of three-jet qqg events vary from rch. = l-10 ± 0.03 (stat. + syst.) [|l^ for the cone jet finder 
to rch.=l-37 ± 0.04 (stat. + syst.) for the JADE-EO jet finder and thus exhibit a strong 
dependence on the jet algorithm employed. 



8 Tests of QCD analytic predictions 

A number of QCD analytic calculations have been presented for the factorial moments of the 
multiphcity distributions of separated gluon and quark jets. Our data permit a test of these 
calculations for the first time. In the following, we test the predictions of analytic calculations 
for the cumulant moments Kg. We study Kg moments, rather than F„ moments, since the results 



for different ranks q are largely independent of each other as was discussed in section 7.2 



An early calculation ^j, valid to the next-to-leading order (n.l.o.) of perturbation theory, 
expresses its results in terms of the strong coupling strength, as, and the number of active 
quark fiavors in the parton shower, n/, allowing the theoretical sensitivity to these quantities 
to be tested. This calculation does not incorporate energy conservation into the parton branch- 
ing processes. More recently, a calculation P] has been presented which is exact for a fixed 
value of as and valid to the next-to-next-to-leading order (n. n.l.o.) if the coupling strength is 
allowed to run. In this paper, we refer to this result as the "n. n.l.o." calculation. By "exact", it 
is meant that the QCD evolution equation is solved without resorting to a perturbative expan- 
sion. Energy conservation, but not momentum conservation, is included in the n. n.l.o. result 
(angular-ordering of partons, introduced to partially account for coherence effects, results in 



approximate momentum conservation [^). The n. n.l.o. results are presented for a fixed value 
of the coupling strength, as=0.22. This value is intended to be an "effective" one, appropriate 
to account for realistic running of as' in the parton evolution of Z° decays. The number of 
active quark flavors is assumed to be nf=4. Therefore, unlike the n.l.o. calculation, the n. n.l.o. 
calculation has not yet been presented in a form which allows the values of as and Uf to be 
varied. The results of the n. n.l.o. calculation are not believed to be strongly dependent on the 



choice of as or Uf or on the use of a flxed as value rather than a running value, however |2^, |28 



In Figures ^ and we present the predictions of the analytic calculations for the Kg moments 
in comparison to our measurements from Figure |^. The results are shown for gluon jets in 
Figure |^ and for quark jets in Figure ^ Besides the n.l.o. and n. n.l.o. results, we show 
the leading order (l.o.) results, obtained from the n.l.o. equations by dropping the n.l.o. 
correction terms. The n.l.o. formulae are evaluated under three conditions: (1) ^/=5 and 



Ag|"^^= 0.209 GeV [||, (2) nf=3 and Ag|"^^ = 0.340 GeV, for which A^"^^ is derived 

from A^ using the prescription relating A^~'^^ to A^"^"* given in ||3^, and (3) n/=5 and 

A^=^^ = 0.209 GeV, with the energy scale at which as is evaluated reduced from Ec.m. to 
-Ec.m./4.0 We take the midpoint between the extreme values found using these three conditions 
as the central n.l.o. result, and deflne a theoretical uncertainty by taking the difference between 
the central and extreme values: the extreme values are in all cases given by condition (1), 
which yields the maximum predicted values of Kg at n.l.o., and (2), which yields the minimum 

^This choice of energy scales is taken from Isl]. 
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predicted values. 



From Figures H and it is seen that the predictions of the l.o. calculation (the star symbols) 
are always well in excess of the data. It is seen that large negative corrections are introduced at 
n.l.o. (the asterisk symbols with uncertainties in Figures ^ and P), which bring the theory into 
agreement with the data for q=2, but which result in even larger discrepancies between data 
and theory for g=4 and g=5 than are observed at l.o. In contrast, the n. n.l.o. calculation (the 
triangle symbols in Figures ^ and ^ is seen to provide a reasonable qualitative description of 
the gluon and quark jet results for all q values. There remain important numerical discrepancies 
between the n. n.l.o. predictions and our data: the n. n.l.o. results for Kq are 0.14, 0.029, 0.0051 
and -0.00042 for gluon jets and 0.32, 0.18, 0.12 and 0.065 for quark jets [|], for ranks g=2, 3, 4 
and 5, respectively. For q=2 and 3, these results are typically a factor of 3 to 6 larger than the 
experimental results given in Tables and ^ for g=4 and 5, the discrepancies are in some cases 
larger and in some cases smaller than this. Nonetheless, it is apparent from Figures || and |] 
that the n. n.l.o. results represent a striking improvement in the theoretical description of the 
cumulant moment data in comparison to the results provided by the lower order calculations. 
This suggests that higher order corrections and energy conservation are essential to obtain a 
reasonable analytic description of gluon and quark jet multiplicity data, similar to what we 
observed in our study of the mean multiplicity ratio r^h. 0- 

Although the analytic calculations do not provide an accurate quantitative description of the 
gluon and quark jet moments, it can be anticipated that certain factors, such as hadronization 
and the dependence of the predictions on n/ or the energy scale, will be common to the gluon 
and quark jet results. Therefore, in Figure |^, we show the analytic predictions for the ratios rxq, 
defined in section [7.2|. For purposes of comparison, the parton level predictions of Herwig and 



Jetset are shown as well. It is seen that the three analytic results, valid to l.o., n.l.o. and n. n.l.o., 
yield almost identical results for tk^- This agreement, in stark contrast to the very different 
predictions which the calculations provide for the individual gluon and quark jet moments 
(Figures H and ||), suggests that the theoretical uncertainties of tk^ are small, and in particular 
that these ratios are only weakly sensitive to the effect of energy conservation and to the values 
of as and Uf. The theoretical uncertainties evaluated for the n.l.o. results are seen to be 
much larger than the experimental uncertainties in Figures |^ and ^, but much smaller than the 
experimental uncertainties in Figure |^, which supports this conclusion. 

For g=2, the analytic calculations predict that the Kq moment of gluon jets is smaller than 
the corresponding moment of quark jets by a factor of about 2.3 (Figure |^. For higher ranks, 
the difference between the gluon and quark jet moments is predicted to be even larger. These 
results are in good agreement with our measurements, as is seen from Figure ^ 

Comparing the results of Herwig and Jetset at the parton and hadron levels (i.e. com- 
paring the open circle symbols to the dashed horizontal lines and the square symbols to the 
solid horizontal lines in Figure it is seen that the hadronization corrections predicted by 
the models are not negligible, especially for q=2 and g=3. The interpretation of the parton 
level Monte Carlo predictions is somewhat unclear, however, for two reasons. First, the Monte 
Carlo simulations implement cutoffs, denoted Qo, to truncate the parton shower at small parton 
invariant masses. For Herwig and Jetset, Qo ~ 1 GeV. We find the parton level Monte Carlo 
results for to be sensitive to Qq. In contrast, analytic predictions for tk^ do not depend on 
a cutoff and in this sense are more reliable theoretically. Second, the mean numbers of partons 
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present at the end of the perturbative shower are small for the Monte Carlo results shown in 
Figure ^ an average of only 4.6 and 2.9 partons are present, respectively, for 41.8 GeV gluon 
and uds quark jet hemispheres generated using Herwig with our tuned parameter set. The 
corresponding results for Jetset are 4.1 and 2.9. These small numbers make it questionable 
whether the parton level Monte Carlo predictions for have much numerical meaning. To 
further investigate this question, we used Herwig to generate 5 TeV gg and uds qq hemisphere 
jets: the resulting mean number of partons present at the end of the perturbative shower was 
56.7 for gluon jets and 30.2 for quark jets. The 5 TeV Monte Carlo jets avoid the two problems 
mentioned above for 41.8 GeV Monte Carlo jets: the condition i^jet >>Qo effectively elimi- 
nates the dependence of the parton level predictions on Qq, while the large parton multiplicities 
make calculation of higher factorial and cumulant moments numerically sensible. The ratios 
determined using the simulated 5 TeV jets were found to be in good agreement with the 
analytic results shown in Figure |^ at both the parton and hadron levels. We conclude that the 
hadronization corrections predicted by the simulations for 41.8 GeV jets are numerically ques- 
tionable and that the agreement between the data and analytic calculations shown in Figure ^ 
is probably not coincidental. 

The OPAL results are based on event hemispheres. The n.n.l.o. results are based on the full 
event multiplicity distributions (in contrast, the definitions employed for the n.l.o. calculation 
correspond to hemispheres ||35|). To assess the effect of this difference, we used Herwig and 
Jetset to evaluate the Kq moments of the full event charged particle multiplicity distributions. 
The full event results were compared to the corresponding hemisphere results shown in Figures |] 
and 1^. Use of full events rather than hemispheres was found to reduce the magnitude of the 
gluon and quark jet Kq moments, typically by 30-60%, further increasing the quantitative 
discrepancy between the data and n.n.l.o. calculation. This difference does not affect the 
qualitative agreement between the data and n.n.l.o. calculation exhibited in Figures ^ and ^ 
since the visible positions of the data are already near zero on the scales of those figures. The 
effect of the difference between hemispheres and full events was found to be negligible for the 
ratios tk^ shown in Figure |^. 



9 Summary and conclusion 



In this study, we have presented measurements of the multiplicity distributions of gluon and 
quark jets. We have determined their mean, dispersion, skew and curtosis values, and factorial 
and cumulant moments. The gluon and quark jets are defined by inclusive sums over the 
particles in gmci. and uds event hemispheres, respectively, with the gmci. gluon jet opposite to 
a hemisphere containing two identified quark jets in e"'"e~ annihilations (the quark jets for the 
ginci. identification are defined using a jet finding algorithm). These inclusive definitions are 
in close correspondence to the definition of jets used for QCD analytic calculations, allowing a 
meaningful comparison of data with theory. Our results for the gluon jet properties are almost 
independent of the choice of the jet finding algorithm, in contrast to other studies of high energy 
(i?jet> 5 GeV) gluon jets. The energy of the jets in our study is about 42 GeV. 

We find the mean multiplicity values of gluon and quark jets to differ by about 50%, in 
agreement with our earlier result but with a substantially reduced uncertainty. We also 
observe differences between the gluon and quark jet skew and curtosis values: the multiplicity 
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distribution of quark jets is observed to be about twice as skewed (asymmetric) as the mul- 
tiplicity distribution of gluon jets, while quark jets are found to have a larger curtosis value 
(curtosis measures the deviation of a distribution from a gaussian shape) than gluon jets. The 
dispersions of gluon and quark jets are found to be the same to within the experimental un- 
certainties. These results are well reproduced by the predictions of QCD parton shower event 
generators. 

We analyze the gluon and quark jet distributions to determine their normalized factorial and 
cumulant factorial moments. These measurements are used to perform the first test of QCD 
analytic predictions of these moments for separated gluon and quark jets. We base our test of 
the anal5d;ic results on the cumulant factorial moments, Kg, since we observe that Kg moments 
of different rank q are largely uncorrelated with each other, unlike the factorial moments, Fq. A 
recent next-to-next-to-leading order calculation which includes energy conservation is found 
to provide a a striking improvement in the theoretical description of the individual gluon and 
quark jet Kg values, in comparison to the description provided by the leading and next-to- 
leading order calculations. 

Our analysis of the Kg moments reveals large differences between gluon and quark jets. For 
rank g=2, the ratio of the cumulant moments of gluon to quark jets is found to be ri^2=0.30 ± 
0.11 (stat.) ± 0.13 (syst.). For rank g=3, ri^3=0.04 ± 0.15 (stat.) ± 0.18 (syst.). The analytic 
predictions for are found to be in quantitative agreement with the data. 
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Table 1: Charged particle multiplicity distributions, expressed in per cent (%), of 41.8 GeV 
ginci. gluon jets and 45.6 GeV uds quark jets. The first uncertainty is statistical and the second 
is systematic. The statistical and systematic uncertainties are correlated between bins. 
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Table 2: The Fg, Kg and Hg factorial moments of the charged particle multiplicity distribution 
of 41.8 GeV ginci. gluon jets. The first uncertainty is statistical and the second is systematic. 
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Table 3: The Fg, Kg and Hg factorial moments of the charged particle multiplicity distribution 
of 45.6 GeV uds quark jets. The first uncertainty is statistical and the second is systematic. 
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Table 4: Correlation matrix between the mean (nch.), dispersion D, skew 7, curtosis c, and 
the Fg and Kg factorial moments of the charged particle multiplicity distribution, for 41.8 GeV 
ginci. gluon jets obtained using the Jetset Monte Carlo at the generator level. Boxes have been 
drawn around correlation coefficients which have magnitudes of 0.90 or larger. 
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Table 5: The ratios rp , and Th of the Fq, Kg and factorial moments of 41.8 GeV gmci. 
gluon and uds quark jets. The first uncertainty is statistical and the second is systematic. 
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0.40 
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Table 6: Differences between the results of the standard analysis and those found by repeating 
the analysis with the systematic changes listed, for the mean charged particle multiplicity value 
(iT-ch.) of 41.8 GeV gincl. gluon jets, for the {rich.) value of 45.6 GcV uds quark jets, and for the 
ratio Teh. between the {rich.) values of 41.8 GeV gmci. and uds quark jets. For item 4, "jet a" 
and "jet b" refer to the two tagged quark jets against which the gmci. jet recoils and "thrust" 
refers to the thrust axis. For items 3, 4, 5 and 10, the larger of the two listed differences was 
assigned as the systematic uncertainty; the smaller of the two differences is given in parentheses 
for purposes of information. 



23 



1 — I — I — I — I — I — I — I — r 



Herwig 5.9, £^^^=41.2 GeV 

^ ^incl. 

— gg event hemispheres 

- - uds event hemispheres 



3 




Figure 1: The prediction of the Herwig parton shower event generator for the charged particle 
multiphcity distribution of ginci. gluon jets from e+e^ annihilations, in comparison to the Herwig 
prediction for gg and uds event hemispheres: (a) on a linear scale, and (b) on a logarithmic 
scale. The jet energies are 41.2 GeV, corresponding to a cm. energy of 91.2 GeV for the 
generation of the e'^e"^ qqginci. events. 
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Figure 2: Corrected distributions of charged particle multiplicity for (a) 41.8 GeV ginci. gluon 
jets, and (b) 45.6 GeV uds quark jets. The total uncertainties are shown by vertical lines. The 
statistical uncertainties are indicated by small horizontal bars. (The statistical uncertainties 
are too small to be seen for the uds jets.) The uncertainties are correlated between bins. The 
predictions of the Herwig and Jetset parton shower event generators are also shown. Numerical 
values for these data are given in Table |I]. 
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Figure 3: (a) The mean, {rich.), and (b) the dispersion, D, of the charged particle multipUcity 
distribution of 41.8 GeV ginci. gluon jets and 45.6 GeV uds quark jets. The total uncertainties 
are shown by vertical lines. The statistical uncertainties are indicated by small horizontal bars. 
(The statistical uncertainties are too small to be seen for the uds jets.) The predictions of the 
Herwig and Jetset parton shower event generators are also shown. 
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Figure 4: (a) The skew, 7, and (b) the curtosis, c, of the charged particle multiphcity distri- 
bution of 41.8 GeV ginci. gluon jets and 45.6 GeV uds quark jets. The total uncertainties are 
shown by vertical lines. The statistical uncertainties are indicated by small horizontal bars. 
(The statistical uncertainties are too small to be seen for the uds jets.) The predictions of the 
Herwig and Jetset parton shower event generators are also shown. 
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Figure 5: The normalized factorial moments of the charged particle multiplicity distribution, 
Fg, for 41.8 GeV gmci. gluon jets and 45.6 GeV uds quark jets. The total uncertainties are 
shown by vertical lines. The statistical uncertainties are indicated by small horizontal bars. 
(The statistical uncertainties are too small to be seen for the uds jets.) The predictions of the 
Herwig and Jetset parton shower event generators are also shown. 
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Figure 6: The cumulant factorial moments of the charged particle multiplicity distribution, 
Kq^ for 41.8 GeV gmci. gluon jets and 45.6 GeV uds quark jets. The total uncertainties are 
shown by vertical lines. The statistical uncertainties are indicated by small horizontal bars. 
(The statistical uncertainties are too small to be seen for the uds jets.) The predictions of the 
Herwig and Jetset parton shower event generators are also shown. 
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Figure 7: The ratios r^, of the cumulant factorial moments Kg of 41.8 GeV gluon and quark 
jets, in comparison to the predictions of QCD analytic calculations and the Herwig and Jetset 
parton shower event generators. The total uncertainties of the data are shown by vertical 
lines. The experimental statistical uncertainties are indicated by small horizontal bars. The 
uncertainties evaluated for the n.l.o. analytic calculation are described in the text. 
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Figure 8: Analytic predictions for the cumulant factorial moments, Kg, of gluon jets, in 
comparison to the OPAL measurements. The uncertainties evaluated for the n.l.o. analytic 
calculation are described in the text. 
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Figure 9: Analytic predictions for the cumulant factorial moments, Kq^ of quark jets, in 
comparison to the OPAL measurements. The inset shows an expanded view. The uncertainties 
evaluated for the n.l.o. analytic calculation are described in the text. 
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